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Abstract. We study the problem of conservation of maximal and lower- 
dimensional invariant tori for analytic convex quasi-integrable Hamiltonian 
systems. In the absence of perturbation the lower-dimensional tori are degen- 
erate, in the sense that the normal frequencies vanish, so that the tori are 
neither elliptic nor hyperbolic. We show that if the perturbation parameter is 
small enough, for a large measure subset of any resonant submanifold of the 
action variable space, under some generic non- degeneracy conditions on the 
perturbation function, there are lower-dimensional tori which are conserved. 
They are characterised by rotation vectors satisfying some generalised Bryuno 
conditions involving also the normal frequencies. We also show that, again 
under some generic assumptions on the perturbation, any torus with fixed ro- 
tation vector satisfying the Bryuno condition is conserved for most values of 
the perturbation parameter in an interval small enough around the origin. Ac- 
cording to the sign of the normal frequencies and of the perturbation parameter 
the torus becomes either hyperbolic or elliptic or of mixed type. 



1. Introduction 

It is well known that in quasi-integrable analytic Hamiltonian systems KAM invariant tori 
are conserved under conditions on the rotation vectors milder than the usual Diophantine 
condition originally introduced by Kolmogorov [30]. A more general condition was intro- 
duced by Bryuno in Refs. [6] and [7], and it is nowadays knowns as the Bryuno condition. 
Among the most exhaustive studies in this direction we cite those by Riissmann (for a recent 
review see Ref. [41]). In some related problems, such as Siegel's problem (in the analytic 
framework), one knows that the Bryuno condition is a necessary and sufficient condition for 
the dynamics to be conjugated to the linear one, as the work by Yoccoz has shown [49]. In 
the case of area-preserving maps the same result likely holds, and for the standard map this 
has been explicitly verified. For Siegel's problem, Yoccoz also proved that a deep relation- 
ship exists between the radius of convergence of the linearising function and the so-called 
Bryuno function. An analogous relationship between the radius of convergence of the conju- 
gating function and the Bryuno function has been found for the standard map by combining 
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the results of Davie [12] with those of Berretti and Gentile [2]. We mention also the work 

by Ecalle and Vallct [13], where it is shown that, under the Bryuno condition, all analytic 
resonant vector fields and diffeomorphisms admit an analytic correction which make them 
linearisable (as conjectured by Gallavotti [16] , and proved under the usual Diophantine con- 
dition by Eliasson [15], hence by Gentile and Mastropietro [24] with techniques more similar 
to those we use in the present paper). Note that for such a problem the rotation vector is 
fixed and no value of the perturbation parameter has to be excluded, so that the problem 
rather simplifies, as all the diflaculties related to estimating the measure of the allowed val- 
ues for the parameters disappear. For instance there is not the difficulty of including the 
correction in the original vector field with a different unperturbed rotation vector, as in the 
case of the KAM theorem for isochronous systems (cf. Ref. [1] for a discussion within the 
formalism used here). Extensions of the Bryuno condition to other contexts, such as the 
reducibility of skew-products, has been partially provided recently by Lopes Dias [34]. 

In this paper we consider a problem of lower-dimensional tori similar to that considered 
in Refs. [38] and [41], with the main difference being that the normal frequencies vanish in 
the absence of perturbation. Such a problem has been explicitly considered in a series of 
papers, such as Refs. [45], [26], [8], [17] and [23]. We refer to the latter for an introduction, 
and for a review of the existing results. In particular we start by considering the same class 
of Hamiltonian systems with d > 2 degrees of freedom considered in Refs. [17] and [23], 
originally introduced in Ref. [44], 

W=iA.A+iB.B + £/(a,/3), (1.1) 

where (a. A) G T"^ x W and (/3, B) G x are conjugate action-angle variables, with 
r + s = d, and • denotes the inner product both in W and in M*. The perturbation /(a, /3) 
is assumed to be real analytic, so that, if we write 

f{oc,f3)=J2^""'"Mf^)^ (1-2) 

there exist positive constants Fq, Fi and kq such that |9^/i/(/3)| < q\FoF^e^'^°^'^^ for all 
e Z"^, all /3 e T* and all q G Z+. For /3q a stationary point of /o(/3) we call ai, . . . , Os the 
eigenvalues of the matrix 9^/o(/3o)- The case of maximal tori is recovered by setting r = d. 
The general case of Hamiltonians describing perturbations of any convex systems will be 
briefly discussed in Appendix A2, even if the full discussion is deferred to Ref. [19]. Here we 
prefer to concentrate ourselves to the simpler model (1.1), in order to distinguish between 
the more relevant features of the renormalisation group techniques and the more technical 
intricacies pertaining rather to problems of spectral analysis and matrix algebra. 

For £ = the system described by the Hamiltonian (1.1) is integrable. Any solution of 
the form (a, /3, A, B) = [olq + u:t, /3o, Aq, Bq), with w = Aq having rationally independent 
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components, fills densely a lower- dimensional torus with rotation vector u. We call A and 

B the non-resonant and resonant, respectively, action variables. With a shift of the resonant 
action variables we can always assume Bq = 0. The so-called normal frequencies, that is 
the frequencies describing the dynamics of the (/3, B)-variables, vanish for £ = 0, so that 
the considered unperturbed torus is neither elliptic nor hyperbolic (nor of mixed type). We 
refer to such a situation by saying that one has a degenerate torus. 

The frequency map A '^(A) is a local diffeomorphism (in our case it is trivially the 
identity), so that the condition B = defines an r-dimensional manifold M^. {resonant 
submanifold) , which is determined by the space of the non-resonant action variable A; we 
call the latter the non-resonant action variable space. 

We are interested in two different problems. 

(A) One can fix the perturbation parameter (small enough) and study for which rotation 
vectors some invariant tori are conserved, in the spirit of the KAM theorem for maximal 
tori, and as done in most of the papers on such a subject, as Refs. [35], [36], [31], [32], [14], 
[38], [9], [41], [5], [28], and many others. 

(B) Either one can look at a lowcr-dimcnsional invariant torus with fixed rotation vector, 
and study the dependence of such a torus on the perturbation parameter. For instance this 
has been done, with the techniques used here, in Refs. [17] and [23]. 

The same twofold program has been followed, under the usual Diophantine condition, in 
Refs. [21] and [22] in the study of the quasi-periodic solutions and of the spectrum for a 
class of two-level systems in a strong quasi-periodic external field. 

About problem (A) we find the analogue of Riissmann's [41] and Poschel's [38] result for 
systems with distinct^ normal frequencies of order 1. In addition, our result applies also to 
the case of non-distinct normal frequencies, provided that they are all different from zero 
and of order e, that is provided degeneracy is removed to first order. In particular this 
means that our result does not follow from the works available in literature: in principle 
one could think to perform a canonical transformation which introduces normal frequencies 
of order ^/\e\, while keeping the perturbation to order e (as explicitly done in Refs. [45] 
and [10]), but in this way the normal frequencies are still required to be distinct in order 
to apply Poschel's result, whereas we do not need such a condition. On the other hand the 
case of possibly non-distinct normal frequencies (of order 1) has been dealt with in Refs. [3], 
[4] , [50] , [47] and [48] only under the usual Diophantine condition. 

Theorem 1. Consider the Hamiltonian (1.1). Suppose (3o to be such that dpfo{f3o) = 0, and 
assume that the eigenvalues of the matrix 9^/o(/3o) are all different from zero (that is ai^Q 
for all i = 1, . . . , s). Let A C M"" be any open set of the non-resonant action variable space. 
Then for any S > there are so small enough and a subset A* C A such that if\e\ < eq and 



We prefer using the word 'distinct', instead of 'non-degenerate', just to avoid confusion, as we are calling 
the normal frequencies 'degenerate' when they vanish for e = 0. 
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iij A* the system described by the Hamiltonian (1.1) admits a lower-dimensional torus of 
the form 

' a = i) + a{ip,f3o,u},e), 
/3 = /3o + b(V,/3o,u;,e), 
' A = u;+(c^-9^)a('0,/3o,a;,e), ^^■'^> 
B= (ui-d^) b(i/',^o,w,e), 

with the functions a and b vanishing at s = 0, analytic and periodic in xj}, and the Lebesgue 
measure of the set A \ A* is less than 6. The parameterisation in (1.3) is such that ip = 
■00 + <*'t describes a linear flow on . In the case of maximal tori (r = d) the subset of 
phase space which is filled by invariant tori has complement whose Lebesgue measure is less 
than C5, for some positive constant C. 

Wc can interpret Theorem 1 by saying that, in the presence of perturbations, the reso- 
nant tori are destroyed in general, but some of them survive. They are determined by the 
stationary points of the potential function /o(/3)- Let /3o one of these stationary points: 
under the (generic) non-degeneracy conditions assumed on the eigenvalues of the matrix 
d^fo{/3o), we can say that, in correspondence of such a point /3o, there is a conserved in- 
variant lower-dimensional torus. The latter is either hyperbolic or elliptic or of mixed type 
according to the signs of the eigenvalues and of the perturbation parameter s: it is elliptic if 
the eigenvalues have the same sign as e, hyperbolic if they have opposite sign with respect 
to £, and of mixed type otherwise. The rotation vector of such a torus will be found to 
satisfy some Diophantine conditions involving also the normal frequencies. In particular we 
shall find that any maximal torus with rotation vector a; which is a Bryuno vector in is 
conserved provided that e is small enough (depending on uj). By Bryuno vector in MJ" we 
mean a vector u) e M' such that 



oo 



„ 2' 



log 



inf 

0<|fc'|<2" 



< OO. 



\U) ■ u\ 



(1.4) 



For further properties of Bryuno vectors we refer to Section 2. Note that in Refs. [38] and 
[41] slightly different (but equivalent) conditions are found for the rotation vectors of the 
surviving elliptic tori, which can be expressed in terms of a suitable approximation function, 
first introduced by Riissmann [40]; cf. the quoted references for further details. 

Also concerning problem (B) our result does not exist in literature, and it represents the 
natural extension of Ref. [17] and [23] to the case of more general rotation vectors. In this 
case wc still need the condition for the normal frequencies to be distinct, as in Ref. [23]. 
Such a condition could be weakened by assuming only that degeneracy is removed to some 
finite perturbation order; cf. Ref. [18]. Though, we shall impose on the rotation vector 
only the Bryuno condition (1.4), a condition much weaker than Kolmogorov's Diophantine 
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condition for maximal tori and Mel'nikov condition for elliptic lower- dimensional tori, as 
usually assumed (cf. Refs. [45], [26] and [27]). 

Theorem 2. Let u: be a vector in W satisfying the Diophantine condition (1-4)- Suppose 
(3q to he such that dpfoi/^o) = 0, and assume that the eigenvalues of the matrix 9^/o(/3o) 
all different from zero and pairwise distinct (that is Oi ^ for all i = 1, . . . , s and Oj ^ Uj 
for all 1 < i < j < s). Then there exists eq and a set £ C {—so, so), with a density point at 
the origin, such that for all e ^ £ there is a, lower- dimensional torus for the system described 
by the Hamiltonian (1.1) with rotation vector u), which can be parameterised as (1.3), with 
tj: gT'^ and the functions a and b vanishing at e = 0, analytic and periodic in ip. 

A density point for £ at the origin means that the relative Lebesgue measure of the set 
£ n (— £, e), that is meas(f n (— £, e))/2e, tends to 1 as e — > 0. We shall say also, in such a 
case, that the set £ has large relative (Lebesgue) measure in (— £o,£o)• 
If e > wc can rcqiiirc ^ Oj for i ^ j only for positive eigenvalues. If all eigenvalues 
are positive then the corresponding torus is elliptic; if all eigenvalues are negative then the 
corresponding torus is hyperbolic. In the first case the allowed values of s form a Cantor 
set in [0,£o) with large relative measure, in the second one all values in [0,£o) are allowed. 
The obvious analogue holds for £ < 0. 

For both theorems we shall give the proof in the case in which all eigenvalues of i9^/o(/3o) 
are strictly positive and £ > 0, which is the difficult case. All the other cases can be obtained 
with trivial adaptations of the proof. Note also that the case of maximal tori can be obtained 
as a byproduct by setting r = d in the following. 

Finally we mention that if do not require that degeneracy of the normal frequencies be 
removed to first order then the problem can become much harder. For instance if no con- 
dition at all is imposed on the perturbation only partial results exist, and only for the case 
s = 1 and u) a Diophantine rotation vector; cf. Refs. [8] and [9] (see also Ref. [18]). 

The paper is organised as follows. In Section 2 we introduce the Bryuno vectors, and we 
briefly review some properties of theirs, which will be used in the forthcoming analysis. Then 
Sections 2 and 3 arc devoted to the proof of Theorem 1 and of Theorem 2, respectively. The 
proofs heavenly rely, both for notations and results, on Ref. [23], and we confine ourselves to 
give full details only for the parts which are really different. In particular the more technical 
aspects arc deferred to Appendix Al. However, by assuming the results of Ref. [23], the 
discussion below becomes completely self-contained. Finally in Appendix A2 we briefly 
discuss how the analysis can be adapted to deal with more general convex Hamiltonian 
systems. 

2. The Bryuno condition 

Given w e set a; = min{|a;i|, ]a;2|}/niax{]a;i|, |w2|}. Let {qn}^=o be the denominators of 
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the convergents of w. 
The Bryuno function B{lj) is defined as the solution of the functional equation [49] 

B{lu + 1) = B{u;), 

B{uj) = -logu + uB{l/oj), ifwG(0,l). 

Define 

lOggn+1 

qn 



(2.1) 



n=0 

Then it is easy to show that B{lv) < oo if and only if -D(w) < oo [49]. 
Given a> e M'^ and n e Z_|_ set 



an{u;) = inf |a) • i/| , (2.3) 

0<|l^|<2'' 

and define the generalised Bryuno funcUon as 

oo 

B(u) = Y —\og^--. (2.4) 

Definition 1. We shall call 03^ = G M' : -8(0;) < 00} the set of Bryuno vectors in . 
For any open set Vl C we call !Br(fi) the set of Bryuno vectors in f2. 

The reason for this terminology relies on the following result. 

Lemma 1. For r = 2 one has uj £ 032 if and only if Dioj) < 00. 

Proof. Given u; € assume for notation simplicity < W2 < wi. Call a = ^2/^1, and 
set sjj = Wia>o, with wq = {^,ct), so that < a < 1 and loga„(w) = logwi + loga„(wo)- 
Consider the sequence of convergents {Pn/Qn}^=i for a [42]; one has l/2qn+i < \aqn —Pn\ < 
1/qn+i, and jwo • I'l > \aqn — Pn\ for all \i'2\ < qn+i, hence for all ji/j < 2g„+i. 

For each n > define r„ and s„ such that 2'^"~^ < 2qn < 2""" and r„ + s„ + 1 = r„+i. 
Hence for all r„ < r' < r„ + ,s„ one has ar'{u}a) — \aqn — Pn\, which implies 

1 / logg„+i \ 1 1 log 2 logg„+i 

r—Tn 

SO that, by using that ^^^Q<ln^ ^ '^^^ obtains that there exist two positive constants 
Ci and C2 such that Z)(u>)/4 — Ci < B{u}) < D{oj) + C2, and the assertion follows. ■ 

The sequence {a„(w)}5^i is non-increasing, so that it converges to monotonically as 
n ^ 00. By taking possibly a subsequence, we can always suppose Q!„+i(w) < a„(a>), 
strictly. 

Definition 2. Set V^=U \ {0}, and define 

n{v) = {n e Z+ : 2""^ < \u\ < 2"} = inf {n e N : \u\ < 2"} (2.6) 
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for any v ^V^. 

For all v^%\ one has, by definition, \<jJ-v\> a„(^)(w) and 2"('')-i < \v\ < 2"(''). 
Definition 3. Given a non-increasing sequence {a* }5^o converging to 0, define 



oo oo 

^* = E^i«gr7^ r;^5]<2"^ (2.7) 

/or p e Z+. One /las r*_,_i > F* /or a/Z p e Z+ 5mc/i t/iaf F* is finite. 

Lemma 2. Lei Q c M'" be an open set, and let {a* decreasing sequence converging 

to zero such that one has B* < oo and T* ~ Cq /or some finite constant Ci^. Call rJ(Co) the 
subset of Bryuno vectors in 53^(^1) such that an{(jj) > ot^ for all n> 1. Then the Lebesgue 
measure of the set 0,'^{Co) = fl \ 0(Co) is bounded proportional to Cq. 

Proof. The measure of the set il.^{Co) can be bounded by 



meas(0'^(Co)) < const. ^ ^ 

„=0 2"-i<|i/|<2 



oo 



(2.8) 

< const. a;2"''2-("-i) < const. F^.^ < const. Co, 

SO that the assertion follows. ■ 
The Diophantine vectors, that is the vectors satisfying the usual Diophantine condition 

|u'-H>|^, (2.9) 

for all e and for suitable positive constants Co and r, are a particular case of Bryuno 
vectors, with a„(tL') > 2^"^Co. In such a case in order to have the convergence of the sum 
in (2.8), hence to apply Lemma 2, one must have r > r — 1, which is the condition for the 
set of Diophantine vectors to have full measure. 

The condition F* = Co motivates us to introduce a new sequence {7„(u;)}^i, with 
7„(a>) = C(^^a„(u>), such that, by setting a* = Co7* and defining 

oo oo 



rp(u^) = E^"(^)2"^ r; = ^7:2"^ = 1, 

r° 00 (2-10) 

^M^E^iogrT^' B log-, 
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one has Tr{u)) > and B{u)) < B for all and all w e 0(Co). 

Note that if |u; • i/| < Co7n(a;) then \i>\ > 2". This is easily checked by contradiction: if 
\u\ < then |w • i/| > Co7„(w). 

3. Fixing the perturbation parameter: proof of Theorem 1 

We follow very closely Ref. [23] (and Ref. [22]), by confining ourselves to show where the 
analysis differs. Also notations which are not defined below are meant the same as in 
Ref. [23]. 

The multiscale decomposition is performed as in Ref. [23], by using the C°° non-decreasing 
function defined as 

with the only difference that now x„ for n > is defined as Xnix) = x(/?~^(7^)~^(t'')a^), 
with /? = 1/4. We set also 'i^nix) = 1 — Xn{x) for n > 0. 

We define clusters and self-energy clusters as in Ref. [23], and we introduce the self-energy 
value Vt{x; e, u)) and the tree value Ya\{9) according to Ref. [23], equations (5.8) and (5.11). 
Of course one has dVy/du; = dxVTdujX + ^^^VT^ and dVr/de = deVr- Note that here and 
henceforth, with respect to Ref. [23] , we are making explicit the dependence of all quantities 
on w, as we are interested also in changing u for fixed e. 

In terms of the self-energy values we can define the self-energy matrices 

n 

A4t^"l(x;£,u>) = ^7WW(a;;£,w), 

p=0 

/ n N (3.2) 

M^^\x;e,u:)={^Xp{^^\x;e,u>))\ ^ VT(a;;£,u;), 
\p=o / Tes^„_i 

where denotes the set of renormalised self-energy clusters of degree k and scale [n] 
(renormalised means that they do not contain any other self-energy clusters). Such matrices 
are formally Hermitian (cf. Lemma 2 in Ref. [23]), so that they admit d real eigenvalues, 
which we denote by aJ"' (a;; e, a>), . . ., a|J*' {x; e, w). 
The matrices A^["l(a;;£,w) in (3.2) can be written as 



M 

^7Wg(x;e,u;) 7Wj3"j(x; e, w) 
where the labels a and /? run over {1, . . . , r} and {r -\- 1, . . . ,d}, respectively. 



M-\x-,eM=[' r^"'"'""^ -^-'.M-^^^-M, (3.3) 
V7W^"i(x;e,u;) Mll{x;e,ui) ) 



8 



With respect to Ref. [23], we slightly change the definition of the propagator divisors for 
n > (cf. Definition 6 in Ref. [23]); see also Ref. [19]. We set 



-1/2 



(3.4) 



and define the propagators as 



where the s elf- energies X^^\e, iv) are defined recursively as A.'"' (e, a;) = a'"'(i/a'" ^\e,uj);e, 



(jj) for n > 1, with A^"^ (e, w) = ea^-r for j = r + 1, . . . , d and A^°' (e, a;) = for j = 1, . . . , r. 
Therefore if a line i is on scale [n], with n > 1, such that g^'^\u) ■ i/^; e, a;) ^ one has 

(3.6) 



min 

i<j<d 



> 



mm 

i<j<d 



< j3 < n- 1, 



and P is chosen in such a way to make uninflucntial the small changes of the propagators 
when shifting the lines in order to exploit the cancellations discussed in Ref. [23] , Appendix 
A3 (see the proof of Lemma 5 below). If a line £ is on scale [0] and gl^l (a; • i/^; e, w) 7^ the 
condition (3.6) has to be replaced with mini<j<d |(w ■ uc)'^ - A^"'(e,a;)| > Cq /S'^ (jq)'^ / AVd. 

Wc define the renormalised expansion for h{ip, l3o,u},e) = (a(i/', /3o, e), b(i/>, /3o, e)), 
by setting 



h{^, /3o, u,,s)=Yl ^''"'^ ^•^^^o, s), h^(/3o, e) = h^, 



h„ = {hu,i, hu,d), K,^ = ^ ^ Val(6'), 



(3.7) 



where the set of trees ,^ is defined as in Ref. [23], Definition 5. 
We shall impose the following Diophantine conditions: 



> C'o7^(i/), 



u;.i.±VAr(£,u>)±VAr(£,u;) 



(3.8) 



> Co7j^(i/) 
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for all i,j = 1, . . . ,d, for all e and for all n > 0. We shall refer to conditions (3.8) as to 
the first Mefnikov conditions (first line) and the second Mel'nikov conditions (second line). 

Definition 4 Given Co G M+ and an open set ft cMJ' call n*(Co) C il the set of Bryuno 
vectors in ri(Co) satisfying all the conditions (3.8). 

Hence the vectors u: E J7*(Co) verify the condition 7n('*^) > Jn ^'•^^ n > 0, with 
Co7n = '^n ^^'-^ t;he sequence {a*}5^Q defined as in Lemma 2, and the first and second 
Mel'nikov conditions (3.8). 



Lemma 3. Call iV„(6') the set of lines in A{0) which are on scale [n]. One has 
Nn{e)<K2-''M{e), M{9)= \u„\, 

for a suitable constant K. One can take K = 2. 



(3.9) 



Proof. First of all note that one can have Nn{9) only if M{6) > 2" ^. Indeed if a line i is 
on scale [n] then there exists i £ {!,... ,d} such that 



\u) ■ vA 



(3.10) 



so that nii'e) > n. Then one must have > 2"('''^)-i > 2""^ hence M{e) > \i>i\ > 2""^ 
thence K2-'^M{e) >liiK>2. 

Then one proves the bound Nn{9) < m.ax.{2-'' K M (9) - 1, 0} for all n > 0, by induction 
on the number of vertices of the tree. The only case which requires a different discussion 
with respect to Ref. [23], Appendix A3, is the one in which the root line £ is on scale [n] 
and exits a cluster on scale [ut], which has only one entering line, say on scale [n'], with 
n' > n. In such a case one has ut <n oi course, and, for suitable i and j, 



07 



\o:-u'\-^)^;-'\e,u:) 
where v = vii and v' = v^/, so that, for suitable r],r]' G {±1} 



(3.11) 



a; . («. - I.') + W Af-^' {e, u;) + r,' J A {s, a,) 



(3.12) 



which by the Diophantine conditions (3.8) implies n{i/ — u') > n, hence one finds M(T) > 
|i/ - i/'l > 2"-\ if M(T) = Et)ey(T) I'^bI- Call 6' the tree having £' as root line. Then 
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by the inductive hypothesis 7V„(6') = 1 + JV„(6i') < 1 + max{2-"isrM(0') - 1,0}. If the 
maximum is the bound is trivially satisfied, because in such a case Nn{6) = 1 and wc have 
seen that in order to have a line on scale [n] one needs M{9) > 2"~^. Otherwise one has 
Nn{0) < 1 + 2-''KM{e') - 1 < 2-"ii'M((9) - 1 + (1 - 2-"-KM{T)), where 2-"ii'M(T) > 1 
by the inequality |i/ — > 2"~^, provided that one takes K >2. ■ 

Lemma 4. Call Nn{T) the set of lines in A(T) which are on scale [n], for n < ut- One 
has 

M{T)= |i^o|>2"-^-\ 7V„(r) < i4:2-"M(T), (3.13) 

oev{T) 

with the same constant K as in (3.9). 

Proof. The first bound in (3.13) can be proved by reductio ad absurdum as in Ref. [23], while 
the proof of the second one is based on the same argument used for proving Lemma 3 (cf. 
Ref. [23], Appendix A3, for further details). ■ 

Another difference with respect to Ref. [23] relies in discussing the change of scale of the 
lines when performing the cancellations inside the families J^ti when looking for bounds on 
the entries of the matrices A^'"' (x; e, w). 

Lemma 5. Assume that the propagators g^\x]£,u)) can he uniformly hounded for all < 
p <n—l as 

g^^\x;e,u,)\<K^C^\^;)-''\ (3.14) 

for some p-independent constants K\ and K2. Assume also that e is small enough. Then, 
with the notations (3.3), one has 



M''2l{x;e,uj) 



< Be 



< Be 



-k;i2" 



iin{e^ 



(3.15) 



-Kl2" 2 



for suitable n-independent constants B and ki . 

Proof. Again we only discuss the differences with respect to Ref. [23]. First we show that 
no cancellation is needed for self-energy clusters T with Co7j^(M(T)) — '^1'^ • i^J, if is the 
momentum flowing through the entering line of T. Note that we can extract from the 
self-energy value a factor e-i^oM{T)/i <- g-Ko2"<*'<'^'V8_ jf .^g 2~"logl/Q!* = a„, we 

have lim„_>oo On = (because B* < 00), hence for oi • 1/ small enough e~'^°'^ ' ' < 



(C07: 



(M(T)) 



«;o/»<l„(M 



<c2(7:(M(T))r <i6i^-^i • 



Hence we need the cancellations only for self-energy clusters T with C07, 



n(M(T)) 



> 4\U) ■ u\ 



if u is the momentum flowing through the entering line of T. In such a case one can reason 
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as follows. For any line I G A(T) and for any n < one has, by the Diophantine conditions 
(3.8), llw • I/O] - ^Jxf\e,u:)\ > Co^^^^oy if is defined as 



1*^0 1 



(3.16) 



itiev(T) 

ltJ<0 



On the other hand one has {v^l < M{T), so that C'o7*(t,o) > (^oln{M{T)) > 4|w-i/|. Therefore 
we can bound 



07 



> 



W • V>i\ 



1 

> - 
- 2 



|u>-i.,°|-A/A"'(£,a;) 



(3.17) 



This implies the following. When considering a family J-^t, a line £ G A(r) with momentum 
i>e can be on a scale [n^] such that g^"'^ (w • 1^°; £, w) = 0. But in such a case it is obtained, 
by shifting the external lines of T, from a line with non- vanishing propagator, that is from 
a line for which (3.6) holds. Then, even if the bounds (3.6) can fail to hold, one still obtains 
bounds of the same form with the only difference that is replaced with /3^/4 in the first 
line and with 4/3^ in the second line. In particular for /3 = 1/4 the inequalities (3.10) and 
(3.12) are still satisfied for all self-energy clusters in the family J^t- 

This shows that one can reason in Rcf. [23] to deduce the boimds (3.15), which are of 
algebraic nature, and are due to symmetry properties of the self-energy matrices, that is 
A^[^"l'^(x;e,w) = A^[^"l(-a;;£,u?) and X[^"lt(x; £, w) = A^[^"l(x; £, w). ■ 

Lemma 6. Assume that the propagators g^^l (x; e, w) can be uniformly bounded for all < 
p < n—1 as in (3.14.), for some p-independent constants Ki and K2. Assume also that e is 
small enough. The matrices Ai^-'^^x; e , u:) are differentiable in x in the sense of Whitney, 
and for all x' , x where they are defined one has 

7W[^"](a;';e,a;) - A^I^"l(a;;£,w) - a^A^[^"l(x;£, w) {x' - x) = e"^ o{\x' - x\), 



d,M\^^\x;E,uj] 



<Bs\ 



(3.18) 



for a suitable positive constant B. Moreover for all j = 1, . . . , d and for a suitable constant 
A one has 

a,Af"l(ar;£,w)| < ^£2, (3.19) 
where the derivative is in the sense of Whitney. 

Proof. The proof of (3.18) can be performed as in Ref. [23]; cf. in particular Section 6. 
Then property (3.19) follows from general properties of Hermitian matrices. One can refer to 
Ref. [23], Appendix A4, in the case in which the eigenvalues are all distinct. Otherwise one 
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can apply the results on non-analytic Hermitian matrices discussed in Ref. [29], Chapter 2, 
Section 6: one can rely on Rellich's theorem [39] to deduce differentiability of the eigenvalues 
and on Lidskii's theorem [33] to obtain a bound on the derivative. ■ 

Lemma 7. Assume that the propagators gl^l (x; £, u)) can he uniformly bounded for all < 
p < n — 1 as in (3.14), for some p-independent constants Ki and K2- Assume also that e 
is small enough. The self-energies satisfy for all <p < n and all 1 < j < d the 

closeness property 



< Ab-''^^\^ 



and one has 



)^\x;£,u}) < Amm.{£^ , ex^} , j = l,...,r, 



(3.20) 
(3.21) 



for a suitable positive constant A. 



Proof The proof can be performed as in Ref. [23], by using the bounds (3.15) and the fact 
that the matrices (x; £, w) are Hermitian (sec Ref. [23], Lemma 2). 

Lemma 8. Assume that the propagators g^^l (.x; e, a;) can be uniformly bounded for all < 
p < n — 1 as in (3.14), for some p-independent constants Ki and K2. Assume also that e 
is small enough. If g^"'\x;e,u)) ^ then one has 



mm 
j=i,...,d 



x^ - Xy {x; e, u>) 



1 

> — mm 
2 3=1,.. .,d 



■A["1(£,u;) 



(3.22) 



The same holds i/ gl"! (x; e, u>) = but (3.6) are satisfied with 0^ replaced with /3^/4 in the 
first line and with 4/3^ in the second line. 

Proof. The inequality (3.22) can be proved by induction on n. For n = it is trivially 
satisfied as the matrix AfM {x\ e, a;) does not depend on x. Let us assume that the inequality 
holds for all < p < n — 1, and let us show that then it follows also for p = n. First of all 
note that in this case we can apply Lemma 7, so that one has 



(3.23) 



for some constant Ai (we have used w S QS^ to deduce e < const.Co(7i^)^). By 

hypothesis one has g(l"l(x;£,u>) ^ 0, hence, by (3.6), 

a;2-Af](£,u,)|>Cl(7:)^ (3.24) 

for all 1 < j < (i and some constant C'l. One the other hand we can write for all j' = 1, . . . , 
x^ - X^p^ {x; e, u;) > x^- X^ {e, u;) - X^ {x; e, w) - aJ^' (e, u>) , (3.25) 
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where we recall that, by construction, A^"' (e, w) = A^"' ( y A^" ^' (e, w); e, w). Hence if j' > r 
and a; > we can bound, by Lemma 6 and in particular (3.19), 



A;?1(x;£,u;)-A^1(£,u,) < Ae^ \x\-.J^F%^) 



(3.26) 



where we have used the relations (3.23) and (3.24). If j' > r and a; < we can apply the 
same argument by using the symmetry property that A^?' {—x; e, w) belongs to the spectrum 

if X^p\x;e,u) does (because M^-"'\—x;s,oj) = (A^[-"l(a;;e, w))-^, see Lemma 2, (ii), in 
Ref. [23]; cf. also the comments at the end of the proof of Lemma 5 above). If / < r then 

X^;^ (x; e, w) - X}p^ (s, up) | = | A^?' (x; e, w) | < Asx^ = As (x^ - X^ (e, w)) , (3.27) 

by (3.21). By inserting (3.26) or (3.27) into (3.25), and choosing / as the value of j 
minimising — A^^'(a;;e,a;)|, then the assertion follows. ■ 

Lemma 9. Let u) € f2*(Co) and assume that e is small enough, say \e\ < Eq. Then the 
series (3.7) admits the bound \hu,-y\ < He~'^\'^^\e\ for suitable positive constants k and H. 
One has sq = 0(Co(7mo)^) "^o depending on kq. 

Proof. One can proceed as in Ref. [23]. Here we outline only the differences. As a conse- 
quence of Lemma 8, we can prove inductively that for all n > the propagators with scales 
[n] are bounded proportionally to {ln)~^^^ ^^'^ finds, in particular, K2 = 2. Then 
the product of propagators can be bounded by relying on Lemma 3 for the lines on scale 
[n], with n > rn^, and bounding with {^1^0)^'^'^ tJ^e propagators of all lines on scale [n] for 
n < mo- Therefore for any tree 9 e ©^i/^-y one has 

2fc 00 



n 



7K] 



< (const.) Cq 



■2k 



1 

7* 
mo 



n 



7; 



2Af„(e) 



< (const.)'' Cq-^'^ 



(3.28) 



and one can choose mo = mo(Ko)) so the last exponential is less than e'^oW\/A_ gy ^laking 



use of the bound Ho 



-i^W\3\ < e"'"!"!, this produces an overall factor e"'"''!''!/^. This 
completes the proof, and it gives £0 = 0(Co(7^(,)^). ■ 

Note that for Diophantine vectors satisfying the bound (2.8) one has mo = TO(logl/Ko)) 
and one obtains £0 = 0{Cq), for fixed r and Cq. 
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If we are interested in studying the conservation of a maximal torus with rotation vector 
(jj satisfying the Bryuno condition B{(jj) < oo, wc can use directly the sequence {'yn{^)}^=o 
for the multiscale decomposition, without introducing a further sequence {'^n}n'=Q- Then 
the result stated in Lemma 9 holds with j^noi'*') replacing 7^^. 

An important remark is that, in the case of perturbations which are trigonometric poly- 
nomials of degree N in the bound (3.9) one can bound M{9) < kN, and as consequence the 
product of propagators in (3.28) can be bounded as 



n 



< exp f 2KNk y ^ log 



(3.29) 



which implies £0 = 0(c-4^-B(a;))^ 

can compare this result with the one found in Ref. [25] 
for maximal tori, where a bound of this kind with the factor 4 replaced with the likely optimal 
2 was obtained. With the techniques described in this paper some further work is necessary 
in order to reach the factor 2; cf. for instance Ref. [2]. On the other hand an advantage with 
respect to Ref. [25], which relies on using Lie transforms for Hamiltonian flows, is that our 
techniques apply, essentially unchanged, not only to the case of flows, but also to the case 
of diffeomorphisms, as done for instance in Refs. [2] and [20], where the case of the standard 
map was explicitly treated. 

Lemma 10. For all w, w' e f2*(Co) one has 

M [^"1 {u:' -u-e^uj')- M [^"1 (cj • i/; e, w) - 

d^M^-''-\u3 ■v;e,u3)- {J -u) = e^|i/| o(|w' - u;|), (3.30) 
Sa,Al[-"l(x;£,a;) < Be''\u\, 



for a suitable constant B, and, as a consequence, 



for a suitable constant A. 



< Ae'M 



(3.31) 



The proof is deferred to Appendix Al, which represents the core of the technical part of 
the paper (by assuming the results of Ref. [23] for granted). Note that in fact it is enough 
to prove (3.30), because then property (3.31) follows by general properties of Hermitian 
matrices; cf. the comments in the proof of Lemma 6. 

Lemma 11. The Lebesgue measure of the set fl \ 0*(Co) is bounded proportionally to Cq. 

Proof. Let us start with the flrst conditions in (3.8). We can reason as in Ref. [22], and 
write u} = oivl\v\ + /3, with /3 • fc' = 0. Then we deflne Q!(t), t G [—1, 1], such that F{t) = 
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ait)\u\ ± ^)^^\e,{a{t),p)) = tCoj^^^y so that dF/dt = + 0{^))da/dt = Col*^^^y 
cf. the proofs of Lemma 3 and Lemma 4 in Ref. [22] for further details. Given p > 1 
we define 17 as the sets of a; € satisfying the conditions (3.8) for n < p; we also set 
= For each n, for fixed u and i, we call In{i,i^) the sets of a; € fit""-'^! such that 



either |u> • i/ ± ^^)^r{e,i^)\ < Coj*^,) or |a; • i/ ± VArU^,'^) ± ^J^^i^,^)\ < Coll^^y 
We define in the same way the sets J„(i, u), with the only difference that the width of the 
sets is 2Co7^(^) instead of C'o7*(-j^j. By the closeness property of Lemma 7, there is some 
ni(i/) = 0(loglog l/7*^j,p such that all the sets In{i,v>) fall inside Jni(u) for n > ni{v). 
Therefore for all i/ e Z^, alH = r + 1, . . . , rf, and all n < ni {v>) we have to exclude all values 
of a; e rjl""^] which fall inside the set Jn{i,,u)\ we refer to Ref. [21], Section 7, for details. 
Note that a; €E 55^ implies < Cn(i/), for some constant C. Hence we can bound the 

measure of the set of excluded values by a constant times 

d ni(v) * oo 

const. E ^0 E E ^ < const. sCo^2"('-i)7; log log 1/7: 

i=r-\-l i/^l/' n=l ' ' n=0 

oo 

< const. sCo^n2"(''-i)7;, 

which is bounded proportionally to Co by Lemma 2. 
Analogously one discusses the other conditions in (3.8). Simply one defines F(i) = 

a{t)\i^\ ± \J}lf'{e,{a{t),P))± ^J^^\ejj4t)^ = tC^ll^^y so that again one has dF/dt = 
[1^1(1 + 0{^/e))da/dt = C'o7*(y), and one can proceed as before. ■ 

To complete the proof of Theorem 1 we have to prove the last assertion about maximal 
tori, that is that for r = d most of phase space is filled by invariant tori. 

We summarise what we have found so far. The invariant tori are determined by the 
corresponding rotation vectors w. For w € n,(Co) we can parameterise the invariant torus 
with rotation vector w as 

ra = i/' + a(i/',u;,£), 

\ A = w + (w-a^)a(V',w,e), ' 

with i/> G T*^. Moreover the function a is analytic in (as the Fourier coefficients decay 
exponentially), while it is defined only on a Cantorian set of values us. 

For each value of V we can consider the map w ^ A(u;), defined in (3.33). For e small 
enough in (3.33) one has |A — a;| = |(cj • d^) a('0,u;,e)| < for some constant R. Call 
ri*(Co,(i) the open set obtained from n*(Co) by excluding all vectors within a distance d 
from the boundary of fi, i.e. il^{Co,d) = e r2(Co) : d{(jj,d^l) > d}, with d{u),d^) = 
min^,/gan |a> — and define A*{Cq) as the image of f2*(Co, -R|£|) of such a map (note that 
the latter is not just the inverse of the frequency map, rather it is a perturbation of it). 
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Then the measure of the complement of the action variable space filled by invariant tori is 
given by 



provided that the Jacobian in the last integral is well defined (that is the map a; — > A(u>) 

is smooth enough, at least in the sense of Whitney) and is uniformly bounded. This turns 
out to be the case, as the following result shows. 

Lemma 12. The solutions of the equations of motion h{'ip,j3o,u:,e) are differentiable in uj 

in the sense of Whitney for u) G fi,(Co, -R|e|)- 

Proof. The proof (for any value of r < d, not necessarily r = d) can be performed as 
for Lemma 10, with the only difference that now we have to deal with the renormalised 
expansion for /i^,^^ instead of the matrices M\-^\x]£,(jj), hence with trees instead of self- 
energy clusters. The condition d(a;,9n) > ii|£| yields that the actions variables A remain 
in for all values oi tj}. ■ 

As a consequence we can bound the Jacobian in (3.34) by using (3.33), which gives d^jA. = 
1 + duj{oJ ■ djp) a{il},u3,e), and Lemma 12, which assures that the last derivative (in the 
sense of Whitney) is bounded proportionally to s. Therefore we can bound meas(^^(Co)) 
proportionally to Co by Lemma 11, and by taking Co = 0(-\/|ij) (which is allowed by Lemma 
9), we obtain the last assertion of Theorem 1. Cf. also Refs. [11] and [37], where the usual 
Diophantine conditions were considered in the analytic and differentiable case, respectively. 



In the following we assume that lj is fixed, and that it satisfies the Bryuno Diophantine 
condition B{(jj) < oo, with B{u!) defined in (2.4). Set 7„ = C(7^a„(u;). 

Let £ G (£o/4,eo] and set a|j^' = eag and £00^ = Ao. Define no G N such that Co7„o+i < 
2\/A^<Co7„„. We set 



and, by using the sequence {Tnln^O' proceed as in Section 3, for constructing the multi- 
scale decomposition of the propagators. Though, we define Al"! (a;; £, a>) = AM(a;;£, a>) for 
n < uq. 

The main difference is that we shall need the following Diophantine conditions: 




(3.34) 



4. Fixing the rotation vector: proof of Theorem 2 




(4.1) 




(4.2) 
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for all i,j = 1,. . . ,d, for all i/ € such that n{i/) > no and for all n > uq. We do not 

impose any conditions like (4.2) for n < ng, because for such scales one has |a; • i/j > 
so that we can bound |(a; • u)'^ — A["'(e,u>)| with (u> • 1^)^/2 for alH = 1, d. In the same 
way we have excluded in (4.2) the values of e such that n{i') < no- Hence, at the price 
of adding a factor 2^ in the bound of each propagator, we can confine ourselves to impose 
(4.3) only for i/ such that «,(//) > tiq and for n > riQ. 
Then we can prove the following result. 



Lemma 13. Call Nn{0) the set of lines in A^{6) which are on scale [n]. One has 

Nn{9)<K2-"M{0), M{e)= I'^ol' 

t)ey(e) 

for a suitable constant K . One can take K = 2. 



(4.3) 



Proof. The proof proceeds exactly as for Lemma 3, with the only difference that we have to 
deal in a different way with the lines on scales n < no and those with scales n > no- The 
same was done in Ref. [23]. ■ 

In the same way the following result is proved. 

Lemma 14. Call Nn{T) the set of lines in A{T) which are on scale [n]. One has 

M{T)= |i/t,|>2"^-\ Nn{T) <K2-'^''^M{T), (4.4) 

vev(T) 

for a suitable constant K. 

Then Lemma 5 is replaced with the following one. 

Lemma 15. Assume that the propagators g^\x\£,u)) can be uniformly bounded for all 
< p < n — 1 as (3.14), for some p-independent constant K. Assume also that eq is small 
enough. With the notations (3.3) one has 



< Be 



-Kl2' 



1/2 



iin{£^,£:a;^}, 



(4.5) 



M^^l{x;e,u)) 

for all n and for suitable n-independent constants B and ki . 

Therefore we can prove the following estimates. The proof is the same as for Lemma 6, 
as one easily realizes that it works for fixed values of e and u>. 



Lemma 16. Let (jj satisfy the Diophantine condition (I.4) and assume that s is small 
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enough, say \e\ < sq. Then the series (3.7) admits the bound \hi,^^\ < Ae "^l^lle] for suitable 
positive constants k and A. One has sq = 0(Co(7^(,)^) with mo depending on w and kq. 

With respect to Section 3 the first differences appear when deahng with Whitney exten- 
sions of the matrices M. '-"1 (x; e, w) : indeed now u) is assumed to be fixed, while e is the free 
parameter. We define £no = (eo/4,£o] and for n > no, recursively, f„ = ^n-i \ ^n^ where 
£° is the set of values of e e such that the conditions (4.2) are violated. We define also 
F — n°° F 



Lemma 17. For all n>0 and all s,e' € £„ one has 

A^t^"] {x; e', w) - X'^"! {x; s, w) - SeA^'-"! {x; s, w) (e' - s) 



= eo(\e — £ 



aeAl[-"l(a;;e,a;) 
and, as a consequence, 



(4.6) 



B' < 
B' < 



< B, 



5.(a1"1(£,u;)±A;"1(£,u,)) 



1 < J < 



r+l<j<i<d, 



(4.7) 



for suitable positive constants B and B' . 



Proof. The proof can be performed as for Lemma 10, with the parameter e now playing 
the role of the parameters a;. We do not give the details, which, however, have been 
worked out in Ref. [23]. Again the upper bound (4.7) follows from (4.6); cf. analogous 
comments in the proof of Lemma 6. To obtain the lower bound one has to use also that 
Aj(.t; e, u;) = a^e + 0(£^), with aj ^ and ttj ^ Oj for i, j = r + 1, . . . , rf; again cf. Ref. [23] 
for details. ■ 

Lemma 17 implies that the matrices A^^-"' (x; e, w) can be extended in (0,eo) to smooth 
functions (Whitney extensions). Again a closeness property of the self-energies, which 
reads 



Afl(£,u>)-Af-^1(^,'^) 



< Be 



(4.8) 



follows from Lemma 15. As before the bounds (4.8) can be improved for the first r self- 
energies, and give |A^"'(x;e,a;)| < Aminje^, ea;^}, j = l,...r. What really changes with 
respect to the previous case is the estimate of the set of allowed values of e, which explains 
why we have required the stronger condition on 9^/o(y9o) that its eigenvalues are non- 
degenerate. The following result holds. 

Lemma 18. The. Lebesgue measure of the set (0, £o)\£* is bounded proportionally to some 
value G{eo), with G{e) = o{e). 



Proof. As in the proof of Lemma 12 we start with the first conditions in (4.2). By setting 
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£ = e{t), with t G [-1, 1], and defining F{t) = u! ■ u ± y A^"'(e(f),u;) = tCoj^f^^y one has 



\dF/dt\ = \de^X!-j'ie,u;)\ \de/dt\ = Co-fl^^y so that, by using that (e,u;) = aje + 0{e^), 
one finds |de/df| < BCo\/sj*(^^-^ for some constant B. Again, for fixed i/ and i, by the 
closeness property of the sclf-cncrgies, we can impose only the conditions corresponding to 
the scales up to ni{v) = 0(loglog 1/7*^^-)), at the price of enlarging the sets of excluded 
values (by a factor 2). Hence the measure of the set of excluded values e G f n > no, 
found by imposing the first conditions (4.2) can be bounded by 

d 00 niii/) 

const. J2 ^oV^ 12 J2 '^n 

i=r+l n=no ra'=rto 2''-i<|t/|<2" (4 9) 

00 

< const, s Co Yl 7:2"Moglog — = G{eo), 

where we have used that no is uniquely determined by sq. Therefore we have 

< const. — V j:{u:) 2"'' log log \ 

"oo (4-10) 

which tends to as no — > oo (that is as eq — > 0). 

The estimate of the measure of the set of excluded values s G £^'^~^\ n > no, found by 
imposing the second conditions (4.2) can be obtained by reasoning in the same way. In such 



a case we need a lower bound on ^^(y {e,u})± y A^ (e, u>)), which requires \ai — aj\ > a 
for all i j and for some constant a; cf. Ref. [23], Appendix A2, for a similar discussion. ■ 

Appendix Al. Proof of Lemma 10 

If we consider two rotation vectors w, w' e fl{Co), they are characterised by the respective 
sequences {7n{^)}^=o and {7„(a;')}rT=o- For all n > one has 7„(a>) > 7* and 7„(t<;') > 7*. 
We introduce some shortened notations, by setting 

n 

S„(i/,a;) = []xp(AM('^-'^;e,<^)), 

*„(«., o;) = Xp(AM (u; ■ v- s, u»)) J ^„(AW (u; • u- e, a;)). 
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and 

s-l 



S„,,(i/, u.', 0.) = ( n Xp(A[^l (a;' • e, a;'))) (x. (A^ (u;' • u- e, u;')) - X«(AW ■ v- s, u;))) 

n 

( n x,>(AW('^-i^;£,'^))), 0<s<n, {A1.2) 

p=s+l 
s-l 

= ( n Xp(AW (a;' • I/; £, a;'))) (x.(AW (w' • i.; £, a;')) " Xs{A^'^ (u; • i^; e, u,))) 



p=0 

n-l 

0<s<n-l, 



( H Xp(AW(a;.i.;£,a;)))^„(AW(u,.i,;£,a;)), 

p=s+l 
n-l 

vl/„,„(i., u;', a;) = ( n XpIA'^I (o;' • u; e, u:'))) (^„(A["1 (o;' • i.; e, w')) - M^^''^ ' e, ^))) , 

p=0 

where all products have to be meant as 1 when containing no factor. 
Finally we define the Hermitian matrices 

£i„(iy,a;) = (a;-i/)2-M[^"l(u;-z/;e,a;), {A1.3) 

withX[-"l(a;-i/;£,u;) given by (3.3). In the obvious way one defines also Dn{i', u}'). 
Note that one has 

n-l 
p=0 

for all m e N. 

Lemma Al. One has 

Xn(x') - Xn(x) - h(x) {x' ~ x) + o{\x' - x\) , \h{x)\ < <P {j*J-^\x' - x\, {A1.5) 
for a suitable positive constant 

Proof. We can write 

\Xn{x') - Xn{x)\ < fi-\ll)-^ W - X\ f AtdMt)), (^1-6) 

Jo 

where x{t) = (3''^ {'j*)-'^ {x + t{x' - x)) and \d^x{x{t))\ < const. ■ 

By noting that tpn = I — Xn, we see that Lemma Al yields the same bounds as (A1.6) 
also if we replace Xn with 

Lemma A2. For u},u}' G f2*(Co) assume that the bounds (3.28) hold for all n' < n. Then 
one has 

\\g^^\u:' ■ u- £, u,') - ffH (u, • u- e, u:) - d^g^^^ (u> • u- e, u:) (a;' - a;)|| = (7:)" Vl o{\u:' - u>|), 
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\\d^g^^\u-v-e,u)\\<D{^*J-'\u\, {Al.7) 
for suitable positive constants D and S. 
Proof. By using the definition (3.5) we have 

(w' • i^; e, u') - (u; ■ iv; e, u,) = vl/„(i., u,')D-\u, u') - ^^{u, u})D-\u, w) 

= -*„(!., U>)D-\u, U>') {Dr,{u, U,') - Dr,{u, w)) I?-l(l/, U>) 

n 

+ ^vi/„p(^,u,',a,)D-i(i/,a;'), (A1.8) 
where we can write 

Dn{u,u}')-Dr,{u,u}) (A1.9) 
= (w' • 1/ + u> • iy)((u>' - u>) • I/) - (w' • u- e, w') + [^"1 (w • i^; e, w) 

so that we obtain 

£>„(i.,w')--Dn(i^,w) (ALIO) 
= (2u{u ■ u) - Oa,A^[-"l(w • V, e, w)) • (w' - w) + ((w' - w) • i/)^ + £>! o(|w' - a;|), 

by the assumed estimate (3.28). 
If ^nW^ w') 7^ we can bound the last sum in (A1.8) by 

n n 

*„,p(i^,u'',u,) \\D-\v,u:')\\ < const, (i^)"^ E(7;)"' l^^l 1'^' ' '^l 

< const. (7:)-'|i^||c^'-u;|, (Al.ll) 

where we have used Lemma Al to bound '^n,p{i', w', w), and (A1.4) to perform the sum over 
p = 0, . . . , n. Note that in order to profitably use the bound (A1.5) we have to use that the 
bounds (3.28) and the consequent (3.29) imply analogous boimds also for the propagator 
divisors A^"! (x; e, w), without the factor e^: this follows from the fact that (3.4) defines 
functions which are smooth in e and u). 

Still if ^n{i^, w') ^ we can bound in (A1.8) also the matrices D~^{u, u') and D~^{l>, w) 
both proportionally to (7^)"^, so that (Al.7) follows. 

If '^niy^ = call an and a'^ the eigenvalues with minimum absolute value of Dni^, w) 
and £'„(t',w'), respectively. If |a^| > |a„| we can proceed as in the previous case, and we 
obtain the same bound. 

Finally if \E'„(^',u)') = and \a'^\ < |a„|, we can write 

W (w' • V- e, w') - H (w • u- e, w) = -5^ (u; • u- e, w). (A1.12) 
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Moreover, by the assumed bound (3.28) we have that the difference Dn{i^,u}') — Dn{i^,u:) 
is given as in (Al.lO), so that ||_D„(i>', a;') — Z?„(i/, < 4\v\ \u)' — u>\. Hence the difference 
between the eigenvalues of £'„(i^,u;) and Dn{y,i*>') is bounded by C\v\ — a;|, for some 
constant C; this is again a consequence of Lidskii's theorem. Therefore for |q:„| > 2C\u\ 
w| we can bound |a^| > |a„| — C\v\ \u} — > |a„|/2, and (A1.7) follows once more by 
reasoning as before, whereas for |q!„| < 2C\u\ |u>' — u>| we can bound in (A1.12) 

5W(u...;.,u.) < \\D-\.,u;)\\ < 

I"" I (A1.13) 
< -^2C\i'\ \u3' -u}\< const. (7;)-2|iy| \uj' - oj\, 

and (A1.7) follows also in such a case. ■ 

Now we can prove Lemma 10. The proof is by induction on n. One can reason as in 
Ref. [22] . More precisely one writes 

A^[<"] (w' . I/; e, a;') - TW'^"] (w • u; s, w) (^1.14) 

n 

V p 
p=0 s=0 

where the matrices M^P^{x;e,u)) arc defined as in Ref. [23], formula (5.9). In the first term 
of the sum the difference M^p\u;' ■ v; e, w' can be written as a sum over 

self-energy clusters T of differences of self-energy values Vt(w' • t') — Vt(w • i-*), computed 
with rotation vectors and uj, respectively. The latter can be written as a sum of several 
contributions, each of which is given by the product of a factor A{u}') depending on u)' but 
no on w times a factor^ ^{^) depending on uj but not on w' times a difference of propagators 
Ai{u}', w) = {uj' ■ V] £, u)') - (w • I/; £, w), with 

The difference /S.({Lt}' .uj) can be boimdcd according to Lemma A4, proportionally to 
\v\ |u)' — u;|, by using that ne < n and the inductive hypothesis. 

Moreover the decomposition A{u>')Ae{u}' ,u})B{u:) can be made in such a way that the 
factor A{u)') corresponds to a connected subset To of T, while the factor B{uj) is the product 
of factorising factors Sj(u;) corresponding to subsets Ti of T containing lines preceding the 
lines of Tq (again we refer to Ref. [17] and [22] for details). The only factor which requires 
some care is that corresponding to the subset, say Ti, connected to the entering line of T. 
But this can be easily discussed as in deriving Lemma 4. Indeed one can prove by induction 



^ No relation with the Bryuno function B(u}) introduced in Section 2. 
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(on the number of nodes) that, given a subset Ti with the considered structure, one has 

J^p(Ti) < K 2~P M (Ti) for all p < n (note the absence of the summand —1 with the respect 
the analogous inductive assumption one makes for Nh{0) and Nh{T)). 

Then, by taking into account also the factor (7^^^)"* possibly arising from A({u:',u}), each 
factor Bi{u)) can be bounded by i?^*e~'^l'''l/^, where h is the number of vertices in Ti, i>i 
is the momentum of the line £i connecting Ti to To, and Bi is some positive constant (the 
proof proceeds as for Lemma 9). The factor A(u}') can be bounded in the same way, and, by 
taking into account also the factor (7^^^)"'' possibly arising from Ae{u}',u}) and the factors 
Q-k\ui\/2^ it can be boimdcd by B!^" , where fco is the number of vertices in To, and B2 is 
some other positive constant (again the proof proceeds as for Lemma 9, but with k replaced 
with k/2). By writing A{u}') = A{uj) + {A{uj') — A{uj)) one can iterate the construction 
above for the difference A{u;') — A{u;). The only difference with respect to the previous case 
is that now the factor k/2 is replaced with k/4. 

All the other terms of the double sum in (Al.14) can be discussed in a similar way, by 
relying once more on Lemma A2. We omit the details, which can be worked out as in 
Ref. [22]. 

Therefore the property (3.28) of Lemma 10 is proved. 

Appendix A2. Extensions and generalisations 

The analysis performed in this article applies to more general Hamiltonians of the form 

n = no{I)+ef{I,ip), {A2.1) 

where {I,(p) e A x T"^, with A c M"^, are conjugate action-angle variables, £ is a real 
parameter, and the functions Tio and / are assumed to be analytic in their arguments. We 
assume also convexity on Tio and a non-degeneracy condition on / which will be specified 
later. Here we confine ourselves to sketch the basic arguments: full details will be published 
elsewhere [19]. Note that in principle weaker conditions on the unperturbed Hamiltonian 
could be considered; we refer to Ref. [43] for a survey of results under the usual Diophantine 
conditions. 

We say that a vector u>* e M"^ is s-resonant if it satisfies s resonance conditions • = 0, 
for s linearly independent integer vectors f 1 , . . . , z/g . 

With respect to (A2.1), the Hamiltonian (1.1) has a very special form. Even by considering 
Hamiltonians of the form 

H = il-I + £/(^), {A2.2) 

it is not always possible to reduce to (1.1). In fact if I: R'' is the linear operator which 

transforms a>* into (a>, 0) e M"^ x M*, where w has rationally independent components, then 
the action variables (A, B) are mixed together and also terms of the form AiBj appear. 



24 



Though, it is easy to extend the analysis to such a case. And with a Httle further work, we 
can consider also Hamiltonians with any unperturbed Haniiltonian Ti.^ satisfying a convexity 
property (so that the eigenvalues of the matrix det^iHo are all strictly positive). The 
frequency map I — > w = dHo/dl is a local diffeomorphsism, so that if we fix Iq in such a 
way that the corresponding rotation vector w(I) is s-resonant, we can find an immersed r- 
dimensional manifold M^., with r = d—s, containing Iq, on which the s resonance conditions 
are satisfied. We shall call a resonant manifold 

Under the action of the symplectic transformation given by the lift <S of I, we can pass 
to new coordinates, which wc continue to denote with the same symbols, such that in the 
new coordinates the rotation vector has became (w, 0). For simplicity we still call the 
resonant manifold in the new coordinates. 

As wc are interested in local properties (in the action variables) we can assume that a 
system of coordinates adapted to M^. has been fixed, so that we can write I = (A, B) in 
such a way that B = identifies M^. For e = a motion on is determined by fixing 
B = and A = Ag in such a way that the conjugated angles (a,/3) move according to 
the law (q;,/3) {a + ujt,/3), with uj uniquely determined by Aq. This means that the 
unperturbed lower-dimensional tori can be characterized by the rotation vectors w € R*^ 
depending on the action variables A. Hence for e ^ the Hamiltonian describing the 
system can be written as 



so that the subset of x T*^ whose unperturbed invariant tori can be continued under 
perturbation can be characterised by the set of allowed values of u), provided the map 
A I— > u>(A) is a local diffeomorphism, which is true under our hypotheses. 

When the perturbation depends also on the action variables, as in (A2.3), of course one 
needs both equations for action and angle variables: 



The main difference with respect to the analysis in Sections 3 and 4 is that the propagators 
are of the form (3.4), with (x^ — M^-"^{x; e,uj))~^ replaced with {ix — M^-"'\x;e,uj))~^, 
where we can write M^-"'\x;e,u!) = C + A/'[-"l(a;;£,w), if the matrix C is given by 



H = Ho(A,B) + £/(A,B,a,/3), 



(A2.3) 



' A = -edafiA,B,ct,f3), 
^ B = -edf3f{A.B,a,(3), 

a = aAWo(A, B) + e5A/(A, B, a, /3), 
^ /3 = ^Wo(A, B) + sdsfiA, B, a, /3). 



(A2.4) 



C = 



/ 



dino dAdsHo 

XdAdBHo dlHo 





25 



and the matrix A/"'-"! = A/^I-"'l(a;; £,a;) is such that, by extracting the dominant order, and 
setting a; = 0, one has 

A/'1^"1(0; e, w) = AAM (0; £, w) + 0(£2), 

/ \ 



AAl°l(0;e,w) 




{A2.7) 



-edpdAf -sd^dsf -e&'^f {A2.6) 
ed\J edAdsf edAdpf 

V sdsdAf edlf ed^dpf) 

whereas the other terms depending on x which are not negligible with respect with the 
dominant ones appear as 

/0{e'^x) 0{e^T) 0(e2.x2) 
0{e^x) 

^{e^x) 

V 0(e2,T) 

Some deep relations turn out to exist between the matrices All-"! (x; e, u;) and their 
transposed, if E denotes the standard symplectic matrix. Then, by using these relations, 
one can bound the propagators in terms of the eigenvalues of a suitable symplectic matrix S: 

for the latter, besides d harmless eigenvalues of order 1 (in e and x) there arc r eigenvalues 
proportional to x^, while the other s eigenvalues are of the form x^ — A^''' (x, e, a;) + 0(ex^ + 

0(£2) , with A^''^ (x, £, w) = £aj_r(a;), j = r + 1, . . . , d, if ai(w), . . . , as(a;) are the dominant 
terms of the normal frequencies (which depend also on A, hence on u>, in this case). The 
aforementioned non-degeneracy condition on / is that the functions 0^ (0;) are all strictly 
positive. The dependence on u> does not introduce any further difficulties, and in fact 
Whitney differentiability in w (as it appears in the subsequent iterative steps) would be 
enough to carry on the analysis. 

Hence, the situation is very similar to that which has been considered in the previous 
Sections 3 and 4, and one can reason essentially as there. Of course notations become 
more cumbersome, because one has to keep trace also of the action variables (which cannot 
be any more expressed trivially in terms of the angle variables), and more sophisticated 
diagrammatic rules have to be envisaged; again see Ref. [19] for details. But the basic 
estimates and arguments are the same, and the same conclusions hold. 

Acknowledgments. I am indebted to A. Giuliani for many discussions and comments. 
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